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IN-DNLS, considered here is a countable infinite set of coupled one dimensional nonlinear ordi- 
nary differential difference equations with a tunable nonlinearity parameter, i/. This equation is 
continuous in time and discrete in space with lattice translational invariance and has global gauge 
invariance. When = 0, it reduces to the famous integrable Ablowitz - Ladik (AL) equation. Oth- 
erwise it is nonintegrable. The formation of unstaggered and staggered stationary localized states 
(SLS) in IN-DNLS is studied here using discrete variational method. The appropriate functional is 
derived and its equivalence to the effective Lagrangian is established. From the physical considera- 
tion, the ansatz of SLS is assumed to have the functional form of stationary soliton of AL equation. 
So, the ansatz contains three optimizable parameters, defining width (/3~^), maximum amplitude 
and its position xo)- Four possible situations are considered. An unstaggered SLS can be 

either on-site peaked (a;o = 0.0) or inter-site peaked {xo = 0.5). On the other hand, a staggered 
SLS can be either Sievers-Takeno (ST) like mode {xo = 0.0), or Page(P) like mode (xq — 0.5). It is 
shown here that unstable SLS arises due to incomplete consideration of the problem. In the exact 
calculation, there exists no unstable mode. The width of an unstaggered SLS of either type de- 
creases with increasing u > 0. Furthermore, on-site peaked state is found to be energetically stable. 
These results are explained using the effective mass picture. For the staggered SLS, the existence 
of ST like mode and P like mode is shown to be a fundamental property of a system, described 
by IN-DNLS. Their properties are also investigated. For large width and small amplitude SLS, the 
known asymptotic result for the amplitude is obtained. Further scope and possible extensions of 
this work are discussed. 
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I. INTRODUCTION 

The study of energy localization in nonlinear lattices 
has become an important field of research in nonlinear 
dynamics in the past couple of decades In this con- 
text, the subject of intrinsic localized modes (ILM) has 
drawn a considerable attention as it offers appealing in- 
sights into a variety of problems ranging from the non- 
exponential energy relaxation^ in soHds, to the local 
denaturation of DNA double strandsf^. The subject is 
also an intense field of study in material science, and 
nonlinear optic applications 

The necessary condition for the formation of intrin- 
sic localized modes (ILMs) or excitations in translation- 
ally invariant nonlinear systems is the balance between 
nonlinearity and dispersion. Furthermore, by localized 
it is meant that the amplitude of such modes goes to 
zero at the boundaries of the system, which is taken 
to be infinitely large. In other words, the relevant lo- 
calization length scale is much much smaller than the 
system size length scale. There are two broad classes 
of intrinsic localizations in (1 + 1) dimensional nonlin- 
ear continuous systems 6J. Shape preserving localized 
excitations, arising in nonlinear continuous systems by 
satisfying the above mentioned balancing condition are 
called dynamical solitons[3, |^ |^. Solitons in contin- 
uous nonlinear Schrodinger e quat ion (cNLS) is an ex- 
ample of dynamical solitons [ol llClj|. By solitons we usu- 
ally mean moving shape preserving nonlinear excitations, 
though there can be stationary solitons also. Take for ex- 



ample cNLS. One particular one-soliton solution of this 
equation is a stationary soliton^. Breathers belong 
to the second category of ILM in nonlinear systems [g- 
Breathers are spatially localized time periodic solutions 
of nonlinear equations. They are characterized by inter- 
nal oscillations[i 

I20I I21L |2^ . Again, by breathers we usually imply sta- 
tionary localized excitations in nonlinear systems. How- 
ever, under appropriate conditions, nonlinear systems 
may have moving breathersR|. As for examples, we note 
that breathers can be found in continuous systems, de- 
scribed by sine-Gordon (sG) equation and modified KdV 
(mKdV) equation^. Even in cNLS, the stationary one- 
soliton solution is nothing but a breather^^J- So, the 
distinction between solitons and breathers is not always 
very rigorous. Breathers are, however rare objects in con- 
tinuous nonlinear equations and are usually unstable|l2l|. 
It is important in the present context to note that contin- 
uous nonlinear equations may have Galilean or Lorentz 
invariance. For example, KdV and cNLS are Galilean 
invariant(7|. So, a soliton of some fixed amplitude in 
the cNLS and KdV can be Galileo boosted to any veloc- 
ity. Similarly, sG equation has both stationary and mov- 
ing breather solutions 0, 112 • These two solutions are, 
however, connected by Lorentz transformation^- So, in 
dealing with stationary ILMs, we consider that moving 
frame which is at rest with respect to the ILM. 

However, models describing microscopic phenomenon 
in condensed matter physics are inherently discrete, with 
the lattice spacing between atomic sites being a fun- 
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damental physical parameter. For these systems, an 
accurate microscopic description involves a set of cou- 
pled ordinary differential-difference equations (ODDE). 
Coupled ODDEs are also encountered in the study of 
many important pro blems in optics and other branches 
of science [T^ 123 . 123^ . So, it is pertinent to discuss next 
what features of continuous nonlinear equations are pos- 
sibly destroyed and what novel features can arise from 
the discretization of at least one of the variables, say one 
spatial dimension. 

In the general discrete case, Galilean or Lorentz in- 
variance in relevant dynamical equations may not be 
present at all or may not be transparent at the equa- 
tion level. Consider, for example the ALj^iE^ and the 
N-AL equations |2^. First one is the example of an in- 
tegrable nonlinear differential discrete equation, which is 
often referred to as the integrable discretization of the 
cNLS equation. The other equation provides an example 
of a differential discrete nonintegrable nonlinear equa- 
tion, having solitary wave solutions. Most importantly, 
the existence of solitary waves in the N-AL equation can 
be shown analytically ,2^. The solitary wave solutions of 
these equations have continuous translational symmetry, 
which can be seen from the analytical expression of the 
one-soliton solution of the AL equation. This, in turn 
implies that both the AL and the N-AL equations have 
the Galilean invariance. So, also in case of ODDEs, sta- 
tionarity in the ILM will imply the moving frame which 
is at rest with respect to the ILM. 

The replacement of the spatial derivatives by spatial 
differences in the equation of motion implies the reduc- 
tion of symmetry of the Hamiltonian, for systems exe- 
cuting Hamiltonian dynamics. In general, lowering the 
symmetry means enriching the class of solutions, be- 
cause less restrictions are imposed. Of course, solu- 
tions are also lost by lowering the symmetry - namely 
ones which are generated by higher symmetry Let 
us consider in this context two discrete nonlinear equa- 
tions, the Frcnkel-Kontorova (FK)^22. 21] and the dis- 
crete nonlinear Schrodinger(DNLS) equations p^. These 
are obtained by standard discretization of sG and cNLS 
respectively [22, ,23, 28, 29j. The FK model can be used to 
describe a broad spectrum of nonlinear physically impor- 
tant phenomena, such as propagation of charge-density 
waves, the dynamics of absorbed layer of atoms on crystal 
surfaces, commensurable-incommensurable phase transi- 
tions, domain walls in magnetically ordered structures 
etc|23,|23- On the other hand, to name a few, DNLS has 
been used to model the self-trapping phenomenon in non- 
linear waveguide arrays 23], to investigate a slow coherent 
transport of polarons in (1+1) dimension in condensed 
matter physics 28] , and to study the dynamical phase di- 
agram of dilute Bose-Einstein condensates JJ. We note 
that both of these discrete equations are nonintegrable 
while their continuous versions are integrable. It is rele- 
vant in this context to know that kink and antikink solu- 
tions of sG equation, which is the continuous integrable 
version of FK model are moving topological solitons, and 



they arise due to the balance between nonlinearity and 
constraints originating from topological invariants in the 
systemj^. On the other hand, there exists no steady- 
state solutions for a moving kink in the FK model. What 
we obtain instead is static kinks 22]. To understand this, 
we note that the uniform discretization of space variables 
transforms continuous translational invariance to lattice 
translational invariance. This, in turn leads to a periodic 
arrangement of Peierls-Nabarro (PN) potential [2^ Isof . 
Therefore, while the continuous translational invariance 
leads to zero frequency Goldstone modes in the system, 
discreteness introduces the PN barrier, with the bar- 
rier energy Epm|22|. Due to this potential any moving 
kink radiates phonons and loses energy (Ekink)- When 
^-kink < -E-pN, the kink is trapped in one of the potential 
wells and further loss of energy by the kink by radiation of 
phonons takes it to the bottom of the well. This, in turn 
yields static kinks. Similarly, sG breathers arise due to 
the high symmetry of the equation, and consequently are 
unstable towards perturbation^. As the discreteness 
in space variables act as an external symmetry breaking 
perturbation, even a weak discreteness does not allow os- 
cillating breather modes exist as dynamical eigenmodes 
of the sG chain, and breathers are destroyed by radiation 
of linear waves. In case of DNLS, similar analysis has 
been done in a perturbative frame using AL one-soliton 
solution as the zeroth order approximation p3. Isij . This 
analysis also shows that discreteness introduces a trap- 
ping potential for moving solitons and when discreteness 
exceeds a critical value, solitonic modes are trapped lead- 
ing ultimately to pinned or stationary solitons. 

It is already mentioned that the AL equation is an in- 
tegrable discrete nonlinear equation. More specifically, 
the said equation is a countably infinite set of one di- 
mensional nonlinear ordinary differential difference equa- 
tions. This equation is continuous in time, but discrete 
in space with lattice translational invariance. The ex- 
act one-soliton solution of the AL equation is charac- 
terized by two parameters, namely, (3 € [0, 00) and 
k £ [— TT. 7r]|24L 1251. For each j3, there exists a band 
of velocities, determined by the other parameter, fc, at 
which the soliton can travel without experiencing any 
PN pinning from the lattice discreteness 133]. Consider 
now other nonlinear equations in this series, namely 
the N-AL equ ation p^. the modified Salerno equation 
(MSE)|2l|3a and the IN-DNLS 32]. All these equations 
are nonintegrable extension of the AL equation, contain- 
ing tunable nonlinearities. The N-AL equation is postu- 
lated and investigated to study the effect of dispersive 
imbalance on the maintenance of the moving solitonic 
profile. The importance of this equation lies in its ap- 
pearance in the dynamics of vibrons and excitons in soft 
molecular chains 26, 28]. The solitary wave solutions of 
this equation are also characterized by the same two AL 
parameters, P and k. However, only certain values of k 
are allowed, though (3 can take all possible permissible 
values. At the allowed values of fc, the term which im- 
parts nonintegrability disappears. This, in turn makes 
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the solitary waves transparent to the PN potential, aris- 
ing from the lattice discreteness. For other values of k, 
the initial AL one-soliton profiles are observed numeri- 
cally to leave phonon tails behind, causing both slowing 
down and distortion of the initial profile. Important too 
in this context is an analytical investigation in a per- 
turbative framework of the dynamics of a moving AL 
soliton, described by the N-AL equation. This analysis 
suggests that any moving soliton having energy below the 
PN barrier, induced by the discreteness in the lattice will 
be pinned, yielding thereby stationary solitonsp^. 

The IN-DNLS is a hybrid form of the AL equation 
and the DNLS, again with a tunable nonlinearity, the 
tuning of which switches the equation from the inte- 
grable AL equation to the nonintegrable DNLSfs^- To 
gauge the physical significance of this equation, we men- 
tion the foUowings. This equation is studied to inves- 
tigate the discreteness induced oscillatory instabilities 
of dark solitonsjs^l- Furthermore, a discrete electrical 
lattice where the dynamics of modulated waves can be 
modeled by this equation is studied to investigate the 
modulation instability of plane waves (35i|. In the MSE, 
the usual DNLS is replaced by a modified version of 
DNLS, the ADNLS, which involves acoustic phonons in 
stead of optical phonons in condensed matter physics 
Darlance[28l l3l| . The study of this equation is also im- 
portant in understanding the dynamics of vibrons and 
excitons in soft molecular chains. It is important to note 
that both IN-DNLS and MSE investigate the competition 
between the on-site trapping and the solitonic motion of 
the AL soliton[2^|32,|3g|- So, the dynamics of a moving 
self-localized pulse, like the AL soliton in the framework 
of the IN-DNLS or the MSE wiU be subjected to two im- 
portant effects. First one is the PN pinning arising from 
the lattice discreteness and second one is a nonlinear in- 
teraction potential, trying to trap or detrap the localized 
pulse. The cumulative effect of these two interactions is 
expected to be the collapse of the moving self-localized 
states to stable, but pinned solitons. This has indeed 
been observed in a numerical simulation(32j. From this 
discussion so far, it can be concluded that the sufficient 
condition to see the effect of discreteness on the dynam- 
ics of nonlinear excitations is that the discrete nonlinear 
equations must be nonintegrable. This nonintegrability 
can arise directly from the discretization of the continu- 
ous nonlinear equations or by adding integrability break- 
ing terms to integrable discrete nonlinear equations. 

Two important linear PDEs, which play very im- 
portant roles in physics in linear systems are free- 
particle Schrodinger equation and the wave equation 
respectively [33. These equations are, of course used to 
describe dynamics in continuous systems. The eigenvalue 
spectra of these equations are continuous function of a 
parameter, k, called wave vector, with the lim — 00 
and the lim = . In case of systems, described 
by Schrodinger equation with a single-particle potential, 
an attractive potential will create localized states below 
the spectra and these are called 'bound states' of the 



system!^. Furthermore, in one (1 -I- 1) dimensional 
systems, even an infinitesimally small attractive poten- 
tial will create an exponentially localized bound state. 
On the other hand, wave equation being second order in 
time, even in (I-fl) dimension no attractive potential, 
however large can create bound states. On the contrary, 
one can get resonances from attractive potentials. 

When the continuity in spatial variables is replaced by 
lattice continuity, the continuous spectra of linear PDEs 
fragments into bands. The number of bands will de- 
pend on the number of lattice points in the unit cell. 
When linear substitutional impurities are added to sys- 
tems, described by a discrete Schrodinger equation, spa- 
tially localized states are formed in the gap between 
bands [371 l38l|. We note that for a state to be locahzed 
and stable, it must be in the gap of the spectra. Further- 
more, these states being exact eigenstates of the relevant 
Hamiltonian, are stationary localized states (SLS). For 
finite number linear impurities in (H-f) dimension, it 
can be shown that the number of spatially exponentially 
localized states cannot exceed the number of impurities 
and there must be at least one exponentially localized 
state js^j- On the other hand, almost all states are expo- 
nentially localized in fully disordered (1 -I- I) dimensional 
systems 38]. However, with correlated disorder, it is pos- 
sible to have some delocalized states [s^^^- stead of 
linear impurities, if finite number of nonlinear impurities 
are present, we again obtain SLS in such systems. This 
can be analytically shown in the systems, described by 
the DNLS iM 142]. 

The spatially discrete analog of the continuous wave 
equation is the coupled mass-spring systems, with springs 
obeying the Hooke's law[33, H^]- Here again we get 
bands of eigenmodes, depending on the number of mass- 
spring units in a unit cell. The lowest band is called 
acoustic branch, which describes the collective motion 
of the masses. Other bands give optical phonons [i^. 
In systems containing finite number of mass impurities, 
only light mass impurities will form exponentially local- 
ized states above the acoustic band in (1 -I- 1) dimen- 
sion. Similar result is also obtained with impurity in 
springs, 3/j. Here also almost all states are exponentially 
localized in totally disordered systems, whether the dis- 
order is in the mass or in the spring or in bothj38j. How- 
ever, no states are obtained below the acoustic branch. 
Most importantly, states around zero frequency remain 
delocalized 38] . In this system also, one can have non- 
linear impurities, either in the spring, or in the on-site 
potential or in both. Any such impurity will produce SLS 
in the system^]- We end this discussion by noting that 
both continuous and discrete linear systems cannot sus- 
tain any localized mode without broken continuous and 
lattice translational invariance respectively. 

A uniform discrete nonlinear system will have lattice 
translational invariance. Similar to continuous nonlin- 
ear systems with translational invariance, nonlinearity in 
discrete systems can also generate localized modes by bal- 
ancing the delocalization effect without requiring broken 



4 



periodicity. Such localized self-organization are the ILMs 
of discrete nonlinear systems. It is important to note that 
ILMs of a discrete nonlinear system are the exact eigen- 
modes of the nonlinear Hamiltonian, describing the sys- 
tem. As in continuous systems, ILMs in discrete systems 
can also be divided in two broad categories, solitons and 
breathersjH |ll|ll[ll|ll|ll[ll|ll|li|23|2l|23. 
In this case also the separation line is not always dis- 
tinct. Consider for example the AL equation. The sta- 
tionary one-soliton solutions of this equation are noth- 
ing but breathers 0) 113 • ILMs are predominantly oc- 
curring nonlinear excitations in discrete nonlinear sys- 
tems. To understand this, we note that stable localized 
modes must always be either below the band or in band 
gaps|23|- So, the discreteness in spatial variables can 
provide a favorable mechanism for the formation and the 
stabilization of ILMs in discrete nonlinear systems by in- 
troducing finite bandwidths, and consequently accessible 
band edges. This, in turn increases the probability that 
the energy of a localized self-organization in a discrete 
nonlinear system will lie in the band gap. Again, the 
band width of a perfect linear discrete system depends 
on the magnitude of the inter-site coupling term. In a 
single band model, if such coupling is weak, we have a 
narrow band. A discrete nonlinear system with narrow 
bands is called anti-integrable. Such anti-integrable non- 
linear systems are then expected to sustain nonlinearity 
induced localized modes in the band gaps by the above 
argument. There is indeed a mathematical proof of this 
in the literature|3 ^M- Of course, it is not necessary 
to have anti-integrable systems to have breathers. The 
existence of a breather solution in the N-AL equation 
has been shown [26j. In fact, in contrast to continuous 
nonlinear systems, in any general discrete nonlinear sys- 
tems, particularly in nonintegrable systems, stationary 
breathers are predominantly occurring ILMs. 

When localized states are formed below the lower band 
edge of a band, they are unstaggered or symmetric local- 
ized states. These states are symmetric under reflection 
through the center, and of course low energy localized 
modes of the system[2^ [s^l- Furthermore, these sym- 
metric localized states can have its peak at a lattice site 
or in between two lattice sites. The first one is called on- 
site peaked unstaggered localized modes. The other one 
is called inter-site peaked unstaggered localized modes. 
When localized states are formed above the upper band 
edge of a band, they are staggered or antisymmetric lo- 
calized states p3^. These states are antisymmetric under 
reflection and high energy excitations of a system. In 
case of staggered localized states, we analogously have 
odd-parity Sievers-Takeno mode (ST) as well as even- 
parity Page(P) mode'300. It is also to be noted 
that staggered localized states have no analog in contin- 
uous svstems|3^. Another kind of ILMs, called twisted 
localized modes can be found in nonlinear lattices [IS liof . 
In this category also, we can have unstaggered as well as 
staggered localized modes [i^ When these modes 

are stationary modes of the system, they are called sta- 



tionary localized states (SLS). We emphasize again that 
SLS of any type, if they are true eigenmodes of nonlinear 
systems are also discrete breathers fill fl^ . They may be 
called trivial breathers. 

Though it is possible to have in circumstances station- 
ary ILMs in discrete integrable nonlinear systems, sta- 
tionary ILMs are formed mostly in nonintegrable non- 
linear systems. We discuss here the stationary ILMs 
of SLS type. We know that stationary solitons of AL 
equations are examples of SLS in integrable nonlinear 
equations, and these are also breather solutions of the 
same equation. We should, however not fail to note that 
though these breather solutions are band edge states, 
their widths are undetermined. On the other hand, 
the formation of SLS in discrete nonlinear systems de- 
pends critically on two factors, the inter-site hopping 
term which determines the width of bands in the corre- 
sponding linear systems and the strength of the nonlin- 
earity, which determines the energy of the self organized 
localized formation. If the first term is predominant, the 
nonlinearity can produce at best localized modes near 
band edges. When localized states are formed near band 
edges, they are weakly localized. In other words, they 
have large widths and small amplitudes. Since, the move- 
ment of these localized modes does not require large scale 
rearrangement in the lattice, such localized modes can 
be made to move by applying small perturbing fields. 
As the movement of any unstaggered localized state will 
not require an inversion of orientation in any of sites, 
these states can easily move compared to its staggered 
counterpart under small perturbation. Again, in case 
of unstaggered localized states, inter-site peaked states 
will have larger widths and smaller amplitudes compared 
to its on-site peaked counterparts. So, inter-site peaked 
states can be made mobile easily by a small perturbation. 
In the other extreme where nonlinearity is strong, strong 
localized modes having nonzero amplitudes only at a few 
sites are formed. These are, of course high energy ILMs. 
Odd parity Sievers-Takeno (ST) modes and even parity 
P modes in strongly anharmonic lattices are examples 
of such strongly localized modes. Since these modes are 
formed from the acoustic branch of anharmonic lattices, 
they appear above the band and hence are staggered lo- 
calized states. It is further found that ST modes are 
unstable to an infinitesimal perturbation. However, this 
mode is not destroyed by the perturbation. Instead, any 
perturbation makes it move|49j. On the other hand, P 
mode is stable and does not move by small perturbations. 
The mobility difference of these modes can be understood 
by the PN potential. Because of the distribution of am- 
plitudes, ST modes are formed at the the maximum of 
the PN potential and the P modes at the bottom of this 
potential^]. For the P mode to move then we need 
enough energy to excite this mode above the PN poten- 
tial. Consequently, under a perturbation, not sufficient 
to take it out of the well, this mode will remain immobile. 
On the other hand, ST modes being at the maximum of 
the PN potential, no energy is needed to take it out of 
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the well. So, an infinitesimal perturbation can make it 
mobile. The mobility difference of on-site and inter-site 
peaked unstaggered localized modes to an infinitesimal 
perturbation can be understood by the same argument. 

With this background, I plan to study here the for- 
mation of both unstaggered and staggered stationary lo- 
calized states in systems described by IN-DNLSfs^- To 
this end, I plan to examine the dependence of the am- 
plitude and width of the localized modes and also the 
eigenfrequency of these modes on the nonintegrability 
parameter of the equation. The energy of the local- 
ized modes are calculated from the Hamiltonian. To the 
best of my knowledge, a rudimentary asymptotic analy- 
sis of this problem is done using the lattice Green func- 
tion approach [s^. IsH ]. For the detailed study, I plan 
to use the discrete variational approach 23, 42, 52]. In 
nonlinear dynamics, the standard variational approach 
has been applied to continuous nonlinear equations to 
study problems of nonlinear pulse propagation in opti- 
cal fibers, and to soliton dynamicsin massive Thirring 
model, to mention a few [2^ [s^, 13 113 ■ In the discrete 
variational approach, one directly proceeds to search for 
discrete solutions of the coupled discrete nonlinear evo- 
lution equations in a restricted subspace by imposing a 
suitable ansatz for the solution[23]. A procedure of av- 
eraging over the discrete dimensions leads to either a set 
of coupled ode's or a set of coupled algebraic equations 
or both for the solution parameters. Therefore, this ap- 
proach permits one to reduce the dimension of the prob- 
lem from a set of many coupled equations to generally a 
much smaller set of equations, determined by the number 
of parameters in the ansatz to be determined. Clearly, 
this method is advantageous when the number of nonlin- 
ear equations is very large. This method has been ap- 
plied to DNLS, for example to study problems of beam 
steering in nonlinear waveguide arrays [2^, and also to 
understand the formation and stability of static and dy- 
namical solitons in one dimensional systems and Cayley 
trees |43.l5^ . We note in this context that equations like 
DNLS describe the evolution of canonical coordinates of 
the canonical phase spacc^23, .5^. On the other hand, 
AL, N -AL and IN-DNLS in their generic form describe 
the evolution of noncanonical coordinates in noncanoni- 
cal phase spaces 13 125I' 26, 32, 56]. Since, these equations 
are derivable from Hamiltonians, th e ge ometry of the dy- 
namics is automatically symplectic[56j. The noncanoni- 
cal symplectic structure of the dynamics is manifested in 
the structure of the Poisson brackets H I1II2II3I l36l] . 
It is, however, to be noted that there exists a global 
nonsingular coordinate transformation for these equa- 
tions, which transforms the noncanonical coordinates to 
canonical coordinates |36|. Therefore, these equations can 
also be described by canonical coordinates with canon- 
ical Lagrangian and Poisson brackets, having canoni- 
cal symplectic structure |36l l56j . I shall, however pro- 
ceed with the variational procedure with noncanonical 



coordinates. I note that in Hamiltonian dynamics, the 
structure of the Poisson bracket is incorporated in the 
Lagrangian |36l Is^ . But, my analysis is done with the 
appropriate functional, which is also obtainable from the 
Lagrangian. So, the noncanonical symplectic structure of 
the Poisson bracket does not pose any problem of find- 
ing SLS in IN-DNLS. The other side of this analysis is 
the following. It shows how the effective Lagrangian can 
be derived from the knowledge of the Hamiltonian and 
constants of motion using the analogous variational ap- 
proach of finding eigenvalues in standard Sturm-Liouville 
problems^]. In other words, I shall also show that it is 
possible to set up the variational problem for the deter- 
mination of eigenvalues without the prior knowledge of 
the Lagrangian. Finally, I note that it has been seen 
in continuous nonlinear equations that when the vari- 
ational method is applied to analyze solitary wave dy- 
namics, the solitary wave solutions may show instability 
in some range of variational parameters. On the other 
hand, the correct dynamics may not show at all such in- 
stability. So, the variational method can produce false 
instabilitiesjS^, E^. This consideration also applies to 
discrete nonlinear evolution equations. However, I do 
not encounter any undesired instability in my solutions, 
which can be ascribed to the variational method. So, this 
aspect, even though important is not dealt with here. 

The organization of the paper is as follows. In the for- 
mulation section below we present the basic equations to 
be studied. In the next section we present a set of results, 
coming from one particular formulation. In this section 
we also show that our formulation gives exact stationary 
localized states of the AL equation. In the next section 
we present another alternative formulation of the same 
problem. We then present the corresponding results. Fi- 
nally, we summarize our main results in the summary 
section. Besides, this paper contains three important as 
well as relevant Appendices. 



II. FORMALISM 

A. General derivation of the nonlinear IN-DNLS 
equation and the variational formulation of the 
corresponding eigenvalue problem 

We consider a dynamical system having 2N general- 
ized noncanonical coordinates, (t>n}, n — 1, ...N 
in a symplectic manifold 56] . Let U and V be any two 
general dynamical variables of the system. Any sym- 
plectic manifold has a natural Poisson bracket struc- 
ture, defined in terms of the inverse of the symplec- 
tic structure function's^]. So, we now define the fol- 
lowing noncanonical Poisson bracket to characterize the 
manifoldHUll^. 
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<rr ^ ,dU dV dV dU ^ ,^ ,, 

We now consider the following Hamiltonian, H. 

H = (0* c^n+i + 4>Ui ^n) ~ 2iy + 2j.^ln[l + \^,,\^]; (2.2) 



which is obtained from the original IN-DNLS Hamilto- 
nian, H through the transformations, 0„ — > ^^JI (pm n Cz 
Z and ^ ^ m, 113 ■ The corresponding Lagrangian C 
in the scaled variables Is^ is 

C - \ Y.(^nK ~kK) ^"'^ 1^ '/"''^ - H. (2.3) 

71 I I 

The dynamical evolution of the n-th generalized coor- 
dinate, can then be obtained by using Eq.(2.1) and 
Eq.(2.2). 

* 0n = (1 + 

= -(1 + |0„P)(0„+1 + 0„_l) 

- 2i.|0„|2^„, (2.4) 

for n G zls^H^- The other set of equations is obtained 
by conjugation. The same equation can be obtained from 
the Lagrangian by using the standard Lagrangian equa- 
tions of motion. We note that under the global gauge 
transformation, (/)„ 0„e*", Eqs.(2.2), (2.3) and (2.4) 
remain invariant. It can also be shown from Eq.(2.4) 
that M = X^Ti + l<^nP] is ^ constant of motion|3^. 
We now assume that 0^ = A" '^n exp (— i u; t), n G Z 
where A = ± 1. Furthermore, 'i'n, n G Z are taken 
realjUl. Then from Eq.(2.4), we get 

{n *)„ - LO^r. + ) (*«+l + *„-l) 

-I- 2 *3 ^ 0. (2.5) 

This is a nonlinear eigenvalue problem and its solutions 
give frequencies of stationary localized states of IN-DNLS 
eauation|32j. Introducing the above ansatz for 0„, n € 
Z in A/", and H, we get 

M = Y,\n[l + ^l], (2.6) 

n 

H = -2 A ^«'„ *„+i-2i.^*2 ^ 2 i.AA 

n n 

= Ho -f 2 1/ A/". (2.7) 
We define next 

F = H - AA/" (2.8) 



where A is the Lagrange multiplier [s^- Setting 5F — 0, 
we get back Eq.(2.5), when A = w. It is also important 
to note that the functional, F can also be obtained from 
C after introducing the ansatz. In Appendix A, I plan to 
discuss the importance of the functional, F. 

B. Variational approach with sech ansatz 

We first note that the system described by IN-DNLS 
equation, Eq.(2.4) has lattice translational invariance. 
So, this system can only form ILMs, arising from the 
competition between the localizing nonlinearity and the 
dispersion from the inter-site hopping^. As the corre- 
sponding linear system is a discrete single band system, 
this further enhances the propensity of formation of ILMs 
either below or above the band. According to the the- 
ory of localization, any self-localized state in 1-d systems 
will have exponential localization in the following sense. 
The amplitude, of the localized mode at the n-th 
site will show exponential decay with \n\ for large val- 
ues of |n||37L 138^. We should also keep in mind that a 
modulus function (|... I) cannot appear in physical prob- 
lem in its generic form. This type of functions can only 
be obtained in any physical problem in the asymptotic 
limit. Furthermore, when v = , Eg. (2. 4) becomes 
the well known AL eauation (241 l25l l3^ . The one-soliton 
solution of Ablowitz-Ladik(AL) equation can either be 
static or dynamic. For both cases, it has the sech pro- 
file, which satisfies also the other requirement for local- 
ized states in one dimension. So, we use the ansatz, 
'^n — ^ — nrr r , n g Z. This ansatz has also 

cosh p(n — xq) ' 

been used in the previous analvsis|32j. For on-site peaked 
and ST like localized states, xq = 0, and for inter-site 
peaked and P like states, = ± i [2II4I lillirt . We 
further write = 5". While (3~^ gives the half- width 
of localization, <I> denotes the maximum amplitude of the 
states. Now, introduction of this ansatz in the functional 
F makes it an algebraic function of the parameters of the 
ansatz, 

.F(*,/?,A,xo) = H(«',/3,A,a;o) - AAA(«',/3,xo), (2.9) 

and we need to find relative extrema of F with respect to 
variables, 5" and /^[s^. The finding of relative extrema 
with respect to these two variables, 4' and (3 means that 
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dF = should imply the following equations 5£ 



5* 



(2.10) 
(2.11) 



dl3 dl3 
where Ai = A — 2v. For what follows next, we assume 

dN_ 
0* 



that §1^0. Then from Eqs.(2.10) and (2.11), we find 



that 



8Ho 

K = ui = 2v 



(2.12) 



n — — oo 



cosh I3{n — xq) 



Xo) 



(2.16) 



n — —oo 



where 



r„(^',/3,xo) = ln[l 



cosh P{n — xo) 



To evaluate Si{P,Xq), S2{P,xo) and A/"(^, 0, Xn), we 
make use of the famous Poisson's sum formula^a, Ha 



MSll^ which reads 



J2 f{nf3) ^ -J [1 + 2 Y.'^osi^^)] Jiv). 

(2.17) 

This application yields 

^i(/?'^o) = sd^, (2.18) 
S2{P.xo) = I + I E r.(/3,a:o), (2.19) 

where 



and also 

/(vI/,/3,A,a;o) = {Hq, = 0. (2.13) 

The other required equation is 

A^(*, /?, Xo) = C = Constant. (2.14) 

We note that we have three unknowns, namely 
A, ^ and (3. But, we also have three independent equa- 
tions to solve for these unknowns. Hence, the problem is 
well-posed. 

C. Calculation of Ho and N 

Introducing the expression of nEZ'm Hq and J\f 
we get 

Roi^, f3,X,xo) = -2 SiiP,xo)-2i^^ S2iP,xo), 

(2.15) 

where 

oo ^ 

Si{f3,XQ) = cosh/3(n-xo) cosh (3(n + 1 - xq) ^^ere 



TT S 
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^s{P,xq) = cos27rsxo — — 
smh — 



^*(1 + = ^ arcsinhV* 



27r 



E cos27rsa;o T,{^,(3), (2.20) 



s = 1 



where 



sin [^arc sinh 

r,(M/,/3) = - 



sinh^ 



We now define a function, /i(/3, v, A, xq) 

/l(p, I^, A, Xq) = ^-7— 5 

1+ Ai^ ^ Ss{l3,,y,X,xo) 



, (2.21) 



5'3(/3, i^. A, Xq) = 1 + 2 cos27rsa;o -; — — 

s = 1 

Now, with this definition, we have 



sinh^l^ 



fi{l3,iy,X,xo) 
4 A 



Ho -4 A 

^Hq 

*. . _,,,ipxM. ,,,,, 



/i(/3, J/, A,xo) 
9. 1 



(2.22) 
(2.23) 
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Again from Eq.(2.20), we have 



2 ^ o -A sin^ arcsinhv^) 

A/'(*,/3,Xo) = -5 (arcsinhV*)^ + 4 > cos27rsa;o — 5 -, (2.25) 

f3 ssinhy 



and from Eq.(2.25) we in turn get 

dM 2 



'd^ - ~ ^ (arcsinh\/*)2 



4 TT arc sinh \/* ^ sm ( arc smh \/* ) 

> cos27rsxo 2 

/32 sinh 4^ 



s = 1 ''""^ 

J- arc silili \/^) 7^-5 

s = 1 



4 sill- arc siuh/*) ^2 



— > cos27rsa;o 5 coth— -. (2.26) 

/32 sinh2f« (3 ^ ^ 

The calculation of Eq.(2.25) is given in Appendix B. 



In our variational formulation, in principle Xq is another 
parameter to be determined from the extrema of the func- 
tional, F (Eq.(2.8)). Now, the extremization of F with xq 
inclusive will yield along with Eqs.(2.10) and(2.11), the 
following equation. 



^Hq 

dxo 



Ai 



dxo 



0. 



(2.27) 



But, from Eqs.(2.21), (2.22) and (2.25), it can be easily 
proved that as < |a;o| < 1, xq = 0, ±5. 



III. THE VARIATIONAL FORMULATION 
WITH Af CONSTANT : RESULTS AND 
DISCUSSION 

A. Ablowitz-Ladik limit 

In the Ablowitz-Ladik limit, v = 0. To probe this 
limit, we evaluate relevant functions and their derivatives 
along the curve = sinh^ /?. Along this curve, from 
Eqs.(2.20), (2.25) and (2.26) we have 



Af{^,(3,xo) = 2(3 
dN _ 
W ~ 



sinh/3 cosh/3 



(3.1) 
(3.2) 

(3.3) 



Since, = in this case, we also have from Eqs.(2.21) 
to (2.24) 
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d(3 



_ 4 A 
sinh/3 

= 4 A cosh/3. 



(3.4) 
(3.5) 



We find from Eqs.(3.2) to (3.5) that /(*, /3, A, xq) = 
{Ho, A/'}{/3,*} = 0. Furthermore, from Eqs.(2.12), (3.3) 
and (3.4) we get w = — 2 A cosh/3. The energy, 
E = H = — 4 A sinh/3. Due to positive semi- 
definiteness of M, we get from Eq.(3.1) that j3 should 
also be positive semidefinite. This is consistent with the 
one-soliton solution of Ablowitz-Ladik equation. 

We now consider the case oiv ^ 0. For convenience, 
we define 



ff(/3,aro) = - [1 + 2 



00 



TT S 





, cos27rsa;o , , 

^ sinh — 



(3.6) 



Along the line 5* = sinh^ /3, we find that 



/(*,/3,A,xo) = 
= — 8 u g{P,xo) tanh/3 



{Ho, A/'}{/3_3,} 
din Ao{(},xo) 



(3.7) 



where Ao{/3,xo) = sinh/3 g{/3,xo)- When /3 — > 
0, tanh/3 ^inAoiji.xo) _^ ^, and consequently 
f{^,f3,X,xo) ~ — ^/3^, provided is finite. So, when 
{y 13^) ~ 0(1), \E' = sinh'^/3 is an asymptotic solu- 
tion of a localized state with a large width and a small 
amplitude. Eigenvalue, w and energy, E = H of these 
localized states are 



and 
E 



2v - 2{\ + V Aq) cosh/3 
- 2 A 



(A + 1) /3^ 



(3.8) 



4A/3 - A.p\^^L^ - 1] 



-4A/3- ^/33. 



(3.9) 



So, according to this asymptotic analysis, when z/ ^ 0, 
the nonintegrability parameter, v and the width param- 
eter (3 of the SLS are not independent of each other. 
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B. Stationary localized states from IN-DNLS 



We now consider various mathematical aspects of the 
formation of stationary localized states in IN-DNLS. We 
consider first Eq.(2.14) along with Eq.(2.25). We restrict 
ourselves to /3 > 0, which is necessary to keep positive 
semi-definite. Furthermore, in the following analysis, we 
assume that /? < 1. In this situation, we can ignore 
infinite sums in Eq.(2.21) and in Eq.(2.25). Due to this 
approximation, Eq.(2.25) yields = sinh^ a\/]3 where a 
is a constant, as required by Eq.(2.14). Since, the right 
hand side of Eq.(2.14) is taken to be a number constant. 



dp 



C = 2.0 a"^, we have 
(3. This, in turn gives 



a* 



(3.10) 



Now introducing Eq.(3.10) in Eq.(2.13), we get = 0. 
In other words, permissible values of (3 are determined 
from the extrema of Hq as a function of /3. From the 
functional dependence of Hq , and ^' on [3, we ultimately 
get 



ffi(a,/3) 
52 (a, /5) 
vX 



(3 , ^ tanha-v/^ 

cosh/5 ;= — — 1 



= 1 



sinh/3 

tanha-y/^ 

sinh/3 g2{a,P) 



(3.ii; 



We note that for a given value of the parameter, a, 
[3 is determined by the nonintegrability parameter, v. 
Furthermore, Eq.(3.11) yields two positive values of (3 
as roots, if two conditions, namely v\ > Q and and 
l^^l < |f^ciiticai|- The behavior of the smaller root {(3s) 
as a function of v for A = 1 and a — 0.5 and 0.25 are 
shown in Fig.l In Fig. 2 we present the variation of [3s 
as a function of the parameter, a for various values of 
V >Q. It should be noted from these figures that /3s < 1 
for these values of a and the chosen interval of v. So, the 
neglect of infinite sums in Eqs.(2.21) and (2.25) is jus- 
tified. It is a simple exercise to see from Eq.(3.11) that 
when \v\ ^0, (3s ^ cP' ■ Then, for small values of v the 
asymptotic solution is the AL stationary localized state 
solution. This is a very important result. This asymp- 
totic analysis reveals that this stationary localized state 
solution of IN-DNLS continuously moves to the AL sta- 
tionary localized state solution when v Q from either 
side. It is further important to note from Fig. 2 that for 

we have a ~ ■ I ^ 



a < 1, 
quently. 



1.0 + 1/ A' 



V \ > 0. Conse- 



sinh 



VI. + V \ ' 



But, as for this range 



of argument, sinha; 



we have \/^ 



Vl.O 4- 1/ A ' 

This agrees with the existing asymptotic analysis 32j. 

The variation of the large root, (3i as a function of v\ 
for various values of a is shown in Fig. 3. Again, by 



0.90 



0.60 



irrespective of the value of CQ. 



0.30 



0.00 




-1.00 



FIG. 1: This figure shows the variation of the smaller root, /3s 
ofEq.(3.11) asafunction of the nonintegrability parameter, v. 
Since A = 1, these states are unstaggered stationary localized 
states. Curve I : a = 0.5 and Curve II : a = 0.25. 



comparing Fig.l and Fig. 3, we see that as v increases, 
the large root, /?/, of Eq.(3.11) decreases from oo, while 
the other root, (3s increases from zero. So, for a given 
a, the value of VcrMcai is determined by the inflection 
point of Hg. This then implies that equations to solve 
for Perineal , Sia<^ Veritieai are obtained by setting both 
^ = and = 0. While the first condition gives 

Eq.(3.11), the second condition yields Eq.(3.12) as shown 
below. 



53 (/3) = ^^^y^ ^"^^ 54(/3) 



P 



sinh P ' 



55(a,/3) = 1 - 53(4^2/3) -4 33 (4«'/3)(l- 53(a'/?)), 
ge{a,P) = l-<73(4a' /3) (1 + 4 cosh/3 54(/3)), 
gm 53(«'/3) + 2 5l(/3)), 



97{a,P) 
9&{a,P) 



'02 cosh2aV/3 95{a,P), 



or j — 
gg{a,P) = -2 A -2 cosh2a ^/P gi{P) gQ{a,P), 



a 

P' 
a" 



gio{a,P) = -4: X ^ cosh^ a^/p gj{a,P), 
g&{a,P) + g^ia.P) + gio{a,P) = 0. 



(3.12) 



We again note that A — ±1 and i/ in Eq.(3.12) is given 
by Eq.(3.11). Of course, vX is positive. We find from 
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0.00 0.20 a 0.40 0.60 

FIG. 2: This figure shows the variation of the smaller root, 
/3s of Eq.(3.11) as a function of the parameter, a for various 
values of the nonintegrability parameter, v. Since A = 1, these 
states are unstaggered stationary localized states. Curve I : 
V = 1.05, Curve II : v = 0.75, Curve III : v ^ 0.5, Curve IV 
: V = 0.25, and Curve Y : v = 0.10. Each curve is associated 
with a dotted curve which shows the variation of q:^(1 + uX) 
as a function of a for the corresponding value of u. 



\ 

A 


' 1 ' 1 

- 


\\ 










,1,1 ^ ^ ~^ 



0.00 0.50 1.00 v?i 1.50 2.00 



FIG. 3: This figure shows the variation of the larger root, Pi 
of Eq.(3.11) as a function of i/X. If A = 1, these states are 
then unstaggered stationary localized states. Curve I : a — 
0.1 and Curve II : a = 0.25, and Curve III : a = 0.50. 



into Eq.(2.12) and using the same approximation as used 
for finding the roots of Eq.(3.11), we obtain that 



Eq.(3.12) that when a — > 0, I'criticai — oo. Again, 
when a S> 1, t'criticai ~ 0. The functional dependence 
of Vcrititcai on a is shown in Fig.4. 

The other important case is where vX < 0. This means 
that either we have an unstaggered state with — z/ or a 
staggered state with +1^. In this case if l^X] > 1, both 
roots of Eq.(3.11) are negative. Inasmuch as J^{ip, (3, xq) 
is positive semi-definite, it is easy to see from Eq.(2.25) 
that this is not permissible. For this case, from Eq.(3.11) 
as expectedly we obtain that when i^X — > — 1+, /3 0, 
and when z/A 0—, /3 o? . See both Fig.l and Fig. 5. 
Inasmuch as for a < 1, permissible values of /3s < 1, the 
neglect of infinite sums in Eqs.(2.21) and (2.25) is again 
well justified. The variation of fig as a function of v for 
a = 1.0, 0.75, 0.5, 0.25, and 0.10 are shown in Fig.5. It 

is seen from Fig. 6 that when z^A < 0, a = \J i q^j^ ^ 3 
is a very good approximation 32] . Another important as- 
pect is in Fig. 7, which shows that for a given z/ > 0, the 
staggered SLS (A = —1) has larger width than the cor- 
responding unstaggered SLS (A = 1). So, the SLS for 
z/A < are basically localized states with large widths 
and small amplitudes. As for eigenvalues of these station- 
ary localized states, introducing Eqs.(2.20) and (2.23) 



2v 



2 A 



P sinh2aV73 



sinh P 2aV/3 



(3.13) 



The energy of these stationary states is given by 



E 



H 



-4 A 



4iy [ 



Ho + 2iyAf 
sinh^ a y/j3 
sinh P 
sinh^ a^/P 



a^P 



(3.14) 



P in Eqs.(3.13) and (3.14) is the root of Eq.(3.11). For 
P = Ps and a not too large, these equations are well 
justified. We already noted that when z^ = 0, Ps — a^- 
Furthermore, |z/A| <C 1 and also a ^ 1, Ps ^ o? . We 
obtain the respective limiting results for these cases from 
Eqs.(3.13) and (3.14). 
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15.0 




0.0 0.5 1.0 1.5 2.0 
a 

FIG. 4: This figure shows the variation of I'criticai as a func- 
tion of the parameter, a for unstaggcrcd stationary localized 
states, A = 1. i^criticai is obtained from Eqs.(3.11) and (3.12) 
in the text. Note that the curve separates the two states 
region from the no state region. 




-1.0 v?i -0.5 0.0 

FIG. 5: This figure shows the variation of the smaller root, (is 
of Eq.(3.11) as a function of v\ for u\ <Q for various values 
of the nonintegrability parameter, a. For A = —1, these states 
are staggered stationary localized states. Curve I : a = 1.0, 
Curve II : a = 0.75, Curve III : a = 0.5, Curve IV : a = 0.25, 
and Curve V : a = 0.10. 



C. Stability and position of stationary localized 
states of IN-DNLS 

We now discuss the issue of stability of these station- 
ary localized states. We note first that when u = Q, the 
resulting nonlinear equation is the AL equation, which 
has both unstaggered and staggered stationary localized 
states. These are basically band edge states. Our vari- 
ational calculation correctly produces these states of the 
AL equation, by letting (is o? as \v\ 0. Sec Fig.l 
and Fig. 5. Over and above it suggests another state for 
which Pi = oo. This can be easily seen in Eq.(3.11). It 
is again seen from Fig.l and Fig. 2 that for any a, in- 
troduction of any v, however small, with vX > makes 
[3s > o? . We observe that /3j^ gives the half- width of the 
localized state. So, for these localized states, the half- 
width reduces with increasing v. Since for A = 1 and 
V > Q (v\ > 0) implies that the on-site nonlinear trap- 
ping potential is attractive, any positive enhancement of 
u should reduce the half-width of the SLS by effectively 
reducing the inter-site hopping. Whether a given ly de- 
fines an attractive or a repulsive potential also depends 
upon the value of A. So, the above argument will hold 
good whenever vX > 0. When i>X is positive, unstaggered 
stationary localized states characterized by Ps are stable. 



On the other hand, for = 0, if there is any stationary 
localized state corresponding to /3; = oo, it is a state with 
a peak of infinite height at a given site with a half- width 
of a few sites. Again, we see from Fig. 3 that for any a, 
when i/A increases, /?; decreases. This means that the 
half-width increases. On the other hand, the introduc- 
tion of u with uX > should reduce the half-width, as 
our argument suggests. Hence, stationary localized states 
corresponding to /3; are unstable. These states, if exist 
in this system, will be unstable towards perturbation. 

Consider next the case where vX < 0. In this case, 
we have either staggered localized states for positive v, 
or unstaggered localized states for negative i^. First of 
all, there is only one set of stationary localized states. 
Furthermore, < (is < for — 1 < A < 0. See 
Fig. 5. Since, for i/ > 0, staggered localized states are 
stabilized by increasing the half- width (see Fig. 7), states 
characterized by /j, are stable. For v negative and A = 1, 
or > and A = — 1, the on-site nonlinear potential 
is repulsive. So, the expansion of the half-width with 
decreasing v is energetically favorable (see Fig. 7). 

So far our analysis did not include the effect of xq , the 
position of the peak on the formation of stationary local- 
ized states and their stability. But, this is also an impor- 
tant part of the problem. However, even a semi-rigorous 
investigation of this problem in this formulation requires 
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0.00 



0.25 



0.50 
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1.00 



FIG. 6: This figure shows the variation of the smaller root, 
/3s of Eq.(3.11) as a function of the parameter, a for various 
values of the nonintegrability parameter, v. Since A = —1, 
these states are staggered stationary localized states. Curve 
I ; = 0.1, Curve \\ : v = 0.25, Curve \\\ : v = 0.5, Curve 
W -.v = 0.75, Curve N : v ^ 0.90 and Curve Nl : v ^ 0.95. 
Each curve is also associated as in Fig. 2 with a dotted curve 
which shows the variation of ci^(l + v\) as a function of a for 
the corresponding value of v. 



the analytical solution of 5" as exactly as possible from 
Eq.(2.25). As there is no simple analytical way of solving 
Eq.(2.25) for 4", one can take recourse to approximation 
methods like, the method of successive substitutions^^. 
I shall describe in Appendix C how this method can be 
used to get approximate dependence on xq, of /3, lu and 
E of the SLS. Of course, the other possibility is to find 
real positive roots of Eq.(2.25) graphically. Inasmuch as 
the exact analytical solution of 5* as a function of /3 is 
difficult in this approach, we shall not follow the present 
line of investigation further. On the contrary, we shall 
show next how the exact dependence of the parameter, 
P, the frequency uj and the energy E of the SLS on xq can 
be obtained by a rational alternation in the variational 
procedure. 



IV. THE VARIATIONAL FORMULATION 
WITH Ho CONSTANT : RESULTS AND 
DISCUSSION 

Since, H and A/" (Eq.2.2) and Eq.(2.6) respectively) are 
two constants of motion, from the expression of H, we 
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FIG. 7: This figure compares the variation of the smaller 
root, /3s of Eq.(3.11) for both unstaggered (A = 1) and stag- 
gered (A = -1) stationary localized states as a function of the 
parameter, a for various values of the nonintegrability param- 
eter, V. Curve \ : v\ = - 0.1, Curve II : !/A = 0.1, Curve III : 
v\ = - 0.25, Curve IN ■.v\ = 0.25, Curve V : j/A = - 0.5, and 
Curve VI : i/A = 0.5. 



see that ffo (Eq.(2.7)) is also a constant of motion. So, 
we reformulate in this section our variational problem 
in which, Hq, in lieu of M is taken to be the number 
constant. In this modified variational approach, we take 
F = A2 Hq + 2 A/", where (A2 — 1) is the Lagrange 
multiplier. This modified approach also yields Eq.(2.12) 
for the eigenvalue, w in Eq.(2.5). One other equation 
required to solve for one of the two unknowns, namely 
/3 and is given by Eq.(2.13). These results are also 
derived in Appendix A. From Eq.(2.22) we find that 



(4.1) 



where a is a number constant, yields Hq = — 4Aa, which is 
again a number constant. Most importantly, in this for- 
mulation ^ is determined explicitly in terms of /? within a 
multiplicative number constant, a. Inasmuch as 5* is pos- 
itive semi-definite by definition, the sign of this constant 
should be such that a/i is positive semi-definite. Fur- 
thermore, in this approach Eq.(2.22) yields = 0, 
irrespective of /3. This in turn gives 



8* 



= ^ 



dlnfi{f3,iy,X,xo) 
d(3 



(4.2) 
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where we have used Eqs.(2.23) and (2.24). Eq.(4.2) can 
also be obtained from Eq.(4.1). Now introducing Eq.(4.2) 

in Eq.(2.13), we get -0- = 0. In other words, permissible 

values of f3 are determined from the extrema of A/", as a 
function of /3. The determination of extrema in turn 
needs Eq.(4.2), Eqs.(2.20) and (2.26). 

Before we proceed further, we observe the following. 
Here, we have a variation problem involving two vari- 
ables, ^ and /?. When ^ is expressed as a function of 
(3, we obtain the Hamiltonian, H = H(/?), and SLSs are 
determined from its extrema, which are obtained by set- 
ting ^ = 0. Of course, in stead of (3, we could have 

used ^ as the fundamental variable. Now, when Af is 
constant, the structure of H (Eq.(2.7)) is such that its 
extrema arc determined from the extrema of Hq. We 
have already investigated here this part. On the other 
hand, we also have the option to take Hq constant, as it 
is done in this sections and in sections to follow. In this 
limit the cxtrcma of H arc determined from the cxtrcma 
of Af, provided Af{f}) has extrema. Another equivalent 
way of envisioning the problem comes from Eq.(2.13), 
which is. of course the direct consequence of the struc- 
ture of the Hamiltonian, H (Eq.(2.2)). We can think of 
an effective dynamical system, having two conjugate dy- 
namical variables, /3 and ^. Then, the Poisson bracket 
formula, Eq.(2.13) suggests that the effective or the re- 
duced dynamical system can be described by the Hamil- 
tonian, Ho(4',/3) having a constant of motion, Af{^,j3), 
or vice versa. Stationary localized states in this dynami- 
cal system picture are determined by fixed points (FPs) 
of the effective or the reduced dynamical system. The 
two sets of extrema, obtained from two procedures or 
two pictures may not be identical. So, in the following 
section, 1 investigate this aspect of the problem. 
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FIG. 8: This figure shows two real roots of Eq.(4.3) in the 
text as a function of the nonintegrability parameter, v for 
three values of the parameter, a. Since, A = 1, these states 
are unstaggcrcd stationary localized states. Curves 1(a), 1(b) 
: a — 1.5, Curves 11(a) and 11(b) : a = 1.75, and Curves in(a) 
and in(b) a = 2.0. While (a) or the lower part of all curves 
is for the smaller root, ft, (b) part or the upper part of these 
curves show the variation of the larger root ft . Note that this 
figure also shows the variation of I'crticai as a function of the 
parameter, a. 



A. Equation for the fixed points of the reduced 
dynamical system and results 

Now, if we altogether ignore the infinite sum in 
Eq.(2.25) which defines A''(^', a;o), we obtain 



= sinh [ 



3^7]- (4-3) 



Introducing Eqs.(4.1) and (4.2) in Eq.(4.3), we get the 
equation which determines /3. If we are mostly interested 
in the roots of Eq.(4.3), having magnitude less than imit 
magnitude, we can as before ignore altogether the infi- 
nite sum in Eq.(2.21), which defines /i(/3, z/, A, Xq)- This 
further simplifies the equation, which determines /3. 

When 1/ = 0, using Eqs.(2.21) and (4.1) we find that 
l3s = arcsinha makes the Eq.(4.3) an identity. We shall 
later show that even with the full expression of , 
derivable from Eq.(2.25), the above choice of /3s makes 
identically zero. In other words, irrespective of Af 



or Ho is taken to be a constant in this constrained varia- 
tional approach, we get the same stationary AL solitons 
with u) = — 2Acosh/3 in both cases. Furthermore, when 
v ^ 0(1), we expect from this result that ^ arcsinha. 
This is also borne out in our numerical calculation, albeit 
not shown here. We shall show it in the exact calculation. 

Considering Eq.(4.3), we consider first the unstag- 
gered localized states with > 0. First of all, for ev- 
ery value of a > 0, we find a i^criticah such that for 

> ^'critical, Eq.(4.3) has no real root. On the other hand, 
1^ < ^'critical (a) , wc find two roots of Eq.(4.3) for a given 
value of a > 0. It is also found that Ps is a monotonically 
increasing function of v while /3; is a monotonically de- 
creasing function of v. These features of the solutions are 
shown in Fig. 8. Tlic^u, according to our previous discus- 
sion, stationary localized states characterized by Ps are 
stable, while the states characterized by are unstable. 
In case of staggered localized states having v > 0, we 
find that for \i/X\ > 1, Eq.(4.3) has no root. Further- 
more, < IvX] < 1, Eq.(4.3) has only one root, /3s. We 
find for A = —1, (3s decreases with increasing i/. This is 
shown in Fig. 9. Since, I3~^ gives a measure of the width 
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FIG. 9: This figure shows the variation of the smaUer root, 
/3s of Eq.(4.3) in the text as a function of the nonintegrability 
parameter, v for various values of the parameter, a . Since, 
A = — 1, these states are staggered stationary locahzed states. 
Curve \ : a = 0.5, Curve II : a = 1.0, Curve III : a = 1.5, and 
Curve IV ; a = 2.0. 




FIG. 10: This figure shows the variation of the smaller root, 
/3s of Eq.(4.3) in the text as a function of the parameter, a 
for various values of the nonintegrability parameter, v. Since, 
A = 1, these states are unstaggered stationary localized states. 
Curve I : i/ = 0.75, Curve W : v = 0.40, Curve III : v ^ 0.25, 
Curve IV : v — 0.0. Curve VI is the straight line. Pa = a. 



of the localized states, from our finding we conclude 
that staggered stationary localized states vanish when- 
ever \vX\ > 1. This happens due to effectively repulsive 
on-site nonlinear potential. This potential helps spread 
the amplitude over the whole sample. We have already 
mentioned this. It is also found for both unstaggered and 
staggered cases that for a ~ o(l), a ^ sinh/J^ ~ 13s, and 
hence, in this asymptotic limit, Eq.(4.1) together with 
Eq.(2.21) yields - --^=h==. Fig.lO shows the re- 
sult for the unstaggered states with v > 0. However, the 
result for staggered localized states is not shown here. 
This result agrees with the asymptotic result in Ref.(32). 

The next important aspect is to study the effect of xq 
on the formation of these states. Another equally impor- 
tant aspect is to examine if unstable localized states that 
we find in the truncated equations or equivalently in the 
leading term analysis exist in the exact calculation. We 
first emphasize that the problem can be solved exactly in 
this reformulated version. For our numerical analysis, we 
use the "FindMinimum" program of "MATHEMATICA 
-version 4" . We discuss below the exact solution. 



B. The exact solution 

We consider first the case of unstaggered stationary lo- 
calized states for 1/ > and xq = 0.0. Fig. 11 shows the 
variation of /?root as a function of the parameter, a for 
various values of the nonintegrability parameter, v. The 
corresponding figure to be compared is Fig.lO. We note 
that for V = 0.0, we obtain from Eq.(4.1) and Eq.(2.21), 
0.0, A, ccq) = asinh/?. On the other hand, analyti- 
cally \E'(/3, 0.0, A, xo) — sinh^/3. So, we must have then 
a = sinh/3i-oot- This is clearly obtained in our numerical 
analysis. As = 0.0 for these curves, both the Curve 
I of Fig. 11 and the Curve V of Fig.lO are defined by 
the equation /?root = arcsinha. We further see in Fig. 11 
and also in Fig.lO that for small values of a all curves 
merge simultaneously on the line /3root — o, and Curve I 
(Fig. 11) or Curve V (FiglO). This in turn implies that 
for a ~ o(l), a ~ sinh/3i-oot ~ Pmot- Hence, for on-site 
peaked unstaggered localized states, ^ ^i'o°+ v 
the asymptotic resultjs^- We further note that in the 
exact calculation, we do not find any root corresponding 
to /3; of Eq.(4.3) for any value of a. This conclusion is 
reached from the following observation in our numerical 
analysis. In our numerical analysis, we have used N and 
M number of terms in two infinite sums in Eqs.(2.21) and 
(2.25) respectively. We find that /3; — > oo monotonically 
if both 'N and Af — > cxd, either separately or simultane- 
ously. So, unstable stationary localized states obtained 
from Eq.(4.3) are spurious and due to the truncation er- 
ror. Similar argument should hold good for the analysis 
of Eq.(3.11). 

The variation of /3root of unstaggered localized states as 
a function of v for various values of a from the exact cal- 
culation, and this is shown in Fig. 12. We see that /3root is 
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FIG. 11: This figure shows the variation of /3root as a function 
of the parameter, a, for various values of the nonintegrability 
parameter, v. This is obtained from the exact calculation. 
For this figure. Ho = constant, = 0.0, and A = 1. Curve 
I : i/ = 0.0, Curve W : v ^ 0.10, Curve \\\ : v ^ 0.40, Curve 
VM w = 0.75, Curve N : v = 0.90, and Curve V\ : v = 1.0. 
Curve VII, the dashed curve is the straight line, /3root = a. 



a nionotonically increasing function of v for v > 0.0. We 
note that /^j^^^t gives a measure of the width of the local- 
ized states. We have also argued before that the width of 
the stable unstaggered localized states for this case must 
decrease with increasing v. So, |/3root| must increase with 
increasing v, if it were to characterize stable SLSs. Inas- 
much as /3root satisfies this criterion, stationary localized 
states corresponding to /3root are stable. 

For inter-site peaked unstaggered localized states, hav- 
ing xq = 0.5, the dependence of /3root on a for a fixed v is 
also investigated for various values of ;/ > 0. It is shown 
in Fig. 13. Fig. 13 also includes for comparison the varia- 
tion of /3root as a function of a for xq — 0.0. We note 
that localized states have larger widths for xq = 0.5. 
This, in turn implies that inter-site peaked SLSs will have 
higher energy. Similarly, the variation of /3root as a func- 
tion of v for Xq = 0.5, is shown in Fig. 14. We again 
note that stationary localized states for xq = 0.5 have 
larger widths. Most importantly inter-site peaked SLSs 
show weak dependence on v. 

In Fig. 15 we show the variation of energy of the on- 
site and inter-site peaked stationary localized states as a 
function of i/ for xq = 0.0 and 0.5. When v ^ o(l), the 
solution approximately has the continuous symmetry of 
the solution of = 0.0^J,l23- So, in this limit, both on- 
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FIG. 12; This figure shows the variation of /3root as a function 
of the nonintegrability parameter, v for various values of the 
parameter, a. It is obtained from the exact calculation. For 
this figure. Ho = constant, xq = 0.0, and A = 1. Curve I : a 
= 0.5, Curve II : a = 1.0, Curve HI : a = 1.5, and Curve IV 
: a = 2.0. 



site and inter-site peaked states should have almost the 
same energy. This is clearly seen in Fig. 15. Similarly, on- 
site peaked states are supposed to be more stable than 
inter-site peaked states. This is also clearly seen by com- 
paring Curves 1(b) and 11(b) in this figure. Again, when 
a reduces, the width of the corresponding stationary lo- 
calized state increases. Consequently, the energetic dis- 
tinction between the on-site and inter-site peaked states 
reduces. This is also clearly evident in Fig. 15. 

A comprehensive understanding of these results, delin- 
eating the basic differences of on-site peaked and inter- 
site peaked unstaggered SLSs is definitely required. To 
this end, we note that for ly > 0, these results indicate the 
operation of a nonlinear attractive potential in the sys- 
tem. This effective nonlinear potential is maximally at- 
tractive at lattice sites. From the physical consideration, 
we argue that the attractive potential assumes minimum 
values at the center of any two consecutive lattice sites. 
Since, the system has lattice translational invariance, this 
potential will also have the periodicity of the underlying 
lattice. An attractive potential effectively reduces the 
inter-site hopping of a particle, consequently helping its 
localization. Secondly, any state can be thought of an ef- 
fective particle with an effective mass, executing a motion 
in a potential. From the physical consideration, it is easy 
to see that stronger is an attractive potential, stronger 
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FIG. 13: This figure shows the variation of /Jroot as a function 
of the parameter, a for two positive values of v. This figure 
presents the exact solution for A = 1. For the Curve = 0.0. 
For Curves n(a) and 11(b) v = 0.40, but xo = 0.5 and 0.0 
respectively. For Curves Ill(a) and Ill(b) v = 1.0, but xo = 
0.5 and 0.0 respectively. 



FIG. 14: This figure shows the variation of /3root as a function 
of the nonintegrability parameter, i/ for two positive values of 
the parameter, a. This figure presents the exact solution for 
A = 1. For Curves 1(a) and 1(b) a — 1, but xo = 0.5 and 0.0 
respectively. For Curves 11(a) and 11(b) a = 2.0, but xo = 
0.5 and 0.0 respectively. 



is the localization. Consequently, heavier is the effective 
mass of the particle. Conversely then the inverse of the 
effective mass gives the localization length of the effec- 
tive state. In this picture then the unstaggered SLS for 
1/ > is equivalent to an effective particle sitting either 
at the bottom of any well (a;o = 0.0) or at the top of 
the same well (xq = 0.5). So, the unstaggered SLS with 
xq = 0.0 will corresponds to a heavier effective mass par- 
ticle than the corresponding SLS with xq = 0.5. In terms 
of localization length, the first kind of states will be more 
localized than the second type. Another important de- 
duction from this picture is that energetically, the first 
kind of states should be more stable. These results are 
seen in our numerical analysis. Again, when we increase 
the parameter, a, we increase the maximum amplitude 
of the SLS. In this effective picture, the depth of the 
potential well increases. Similar situation also occurs by 
increasing u. This, in turn implies that the effective mass 
of the particle at the bottom of the well will increase with 
increasing a and v. So, in both cases, the width of the 
on-site peaked unstaggered SLS should decrease, as seen 
in our numerical calculation. The effective periodic po- 
tential, however is a function of at least three variables, 
the position variable, a;o, the parameter, a and the non- 
integrability parameter, v. Since, the SLS with xq = 0.5 
shows weak dependence on i^, our results suggest that 



the top of the potential is not substantially affected by 
the change in v. On the other hand, from our results 
we deduce that tlu^ top of the potential is energetically 
stabilized by increase in a. 

Finally, some exact calculations of Proot for staggered 
SLSs for > 0, A = —1 and for both xo = 0.0 and 0.5 
are presented. In Fig. 16 we present the variation of /3root 
as a function of v for o = 1.0 and 2.0. Consider first 
Xq = 0.0. Curves 1(a) and 11(a) in Fig. 16 are almost 
identical to corresponding curves Fig. 9 with a discernible 
deviation in the magnitude of /3root for large values of a 
together with small values of v. The same calculation 
with Af = constant gives /3s ^ linearly as ^ —1. 
See Fig. 5. Fig. 17 shows the variation of /3root as a func- 
tion of a for various values of > 0. The important 
point to note is that Pmoti'^^ < 0) < l3root{v = 0)- Iri 
other words, SLSs for v > Q are stabilized by expansion 
of the width. As /3root = a fine becomes tangent to all 
curves in Fig. 17, we infer that when a ~ o(l), /3root a. 
Again on the curve z/ = 0, sinh/3root = a. In this asymp- 
totic limit, we then have from Eqs.(2.21) and (4.1) that 
\/^ '^ff=- The corresponding approximate calcula- 
tion for the model with J\f = constant is shown Fig. 6. 
Fig. 6 shows that when a ^ o(l), (3s ^ o?{\ — v). This in 
turn yields the same asymptotic result for Compar- 
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FIG. 15; This figure shows the variation of the energy of un- 
staggered stationary locahzed states as a function of v for 
two values of a and for two permissible values of xq. Of 
course, the result is obtained from the exact calculation with 
Ho = constant. For Curves 1(a) and 11(a) a = 1.0, but 
= 0.0 and 0.5 respectively. For Curves 1(b) and 11(b) 
a — 2.0, but xo = 0.0 and 0.5 respectively. 




FIG. 16: This figure shows the variation of /3root as a function 
of the nonintegrability parameter, u for a — 1.0 and 2.0, as 
obtained from the exact calculation. Since, A — —1, these 
states are staggered stationary localized states. Curve 1(a) : 
a = 1.0, and xo ~ 0.0. Curve 11(a) : a — 2.0, and xo = 0.0. 
Curve 1(b) : a = 1.0, and xq = 0.5. Curve 11(b) : a = 2.0, 
and Xo — 0.5. 



ing our results for staggered SLSs from three different 
approaches, namely, approximate calculations with (a) 
JV = constant, (b) Hq = constant, and (c) the ex- 
act calculation with Hq = constant and xq — 0.0, we 
conclude that all three give qualitatively the same result. 

We now consider the basic difference between stag- 
gered SLS with Xq = 0.0 and staggered SLS with xq = 
0.5. Fig. 16 shows the the variation of /3root as a function 
of 1/ > for a = 1.0 and 2.0. From this figure, we see that 
though /3root ^ as J/ 1 for both xq = 0.0, and 0.5, 
the magnitude of /Sroot for intermediate values of v is 
greater for the SLS with xq — 0.5. This, in turn implies 
that SLS with xq — 0.5 has smaller localization length. 
This is very much opposite to what we observe for an 
unstaggered SLS. We consider next Fig. 18. We define 
AEim = El — Em, where Ei and E^ define the energy 
of staggered SLS with xq = 0.0 and x^ ~ 0.5 respec- 
tively. Fig. 18 shows the variation of AEim as a function 
for ly G [0.1). As Ai?/,„ > for intermediate values of i/, 
SLS with Xq = 0.5 is energetically stable compared to its 
equal counterpart. 

In further analysis, we note that the amplitude distri- 
bution of staggered SLS around the maximum amplitude 
site is ST mode like for xq — 0.0 and P mode like for 



Xq — 0.5. In this case also there is a periodic arrange- 
ment of potential wells in the system. Since, the on-site 
nonlinear potential is repulsive in this case, this potential 
will attain the minimum at the center of two consecutive 
lattice sites. So, the whole periodic arrangement of wells 
is shifted by half a lattice constant. Consequently, the 
effective particles corresponding to the P like mode and 
the ST like mode are sitting at the bottom of a well and 
at the top of a well respectively [H^- So, the P like mode 
corresponds to an effective particle with larger effective 
mass than the corresponding ST like mode, and thereby 
having a smaller localization length. From this picture, 
we also deduce that the P like mode is energetically more 
stable than the ST like mode. Most importantly however, 
we prove by our variational approach that the existence 
of P like mode and ST like mode is a fundamental prop- 
erty of a system described by IN-DNLS. To successfully 
explain the dependence of the localization length of the P 
like modes on these two parameters, a and v respectively, 
we need the following behavior of the potential. When 
the parameter, a increases, the effective well depth must 
increase for intermediate values of v. Consequently, the 
effective mass of the particle will increase, and the local- 
ization length will decrease. But, in case of v, the effec- 
tive well depth must decrease as v I. Furthermore, as 
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FIG. 17: This figure shows the variation of /3root as a func- 
tion of the parameter, a for positive values of v. This figure 
presents the exact solution for A = —1 and xq = 0.0 . Curve 
l:v^ 0.995, Curve II : v = 0.75, Curve III : i/ = 0.40, Curve 
lY : V = 0.1, and Curve Y : v ^ 0.0. Curve VI, the solid 
curve is the straight line, /Jroot = a- 

discussed in the context of unstaggered SLS, the depen- 
dence of the locahzation length of ST hke modes deter- 
mines the dependence of the top of the potential well on 
two important parameters, a and v. 



V. SUMMARY 

IN-DNLS is a one dimensional discrete nonlinear equa- 
tion with a tunable nonintegrability parameter, v 32, 33|. 
When = 0, it reduces to the famous AL equations,245 
|25| . The importance of IN-DNLS in physics as well as in 
nonlinear mathematics is discussed in the text. In this 
paper, primarily eigenvalues, energies and corresponding 
site-amplitudes of SLSs of IN-DNLS are studied using 
discrete variational formulation[23. li^ . l53 |. The stan- 
dard variational approach starts from the respective La- 
grangian to study this type of problem. In this paper, 
however the appropriate functional is derived using the 
standard variational procedure for finding eigenvalues 
of Sturm-Liouville equations 58] . In other words, it is 
shown here how the effective functional can be derived 
from the Hamiltonian and constants of motion without 
the prior knowledge of the Lagrangian. The uniqueness 
of the functional is also established by showing its equiv- 
alence to the effective Lagrangian. 
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FIG. 18: AEim = El — Em, where Ei and Em define the 
energy of staggered SLS with xq = 0.0 and xq = 0.5 respec- 
tively. This figure shows the variation of AEim as a function 
for p e [0.1). Curve I : a = 1.0, and Curve II : a = 2.0. 



Inasmuch as localized states in one dimensional linear 
impure as well as disordered systems show asymptotic 
exponential decay [s^, a "sech"' ansatz with two pa- 
rameters, (3 and <I> = is used to find eigenvalues and 
eigenfunctions. In this choice, (3^^ and define the 
width and the maximum amplitude of SLS respectively. 
This ansatz is so chosen as it gives AL stationary local- 
ized states when z/ ^ 0. Furthermore, SLSs of IN-DNLS 
are assumed to belong to the class of breathers with a sin- 
gle frequency 0. Since, stationary solitons of AL equa- 
tions are breathers of this class, this choice of form for 
SLSs of IN-DNLS is justified. 

Very naturally two procedures have emanated in our 
variational calculation. In the first case, the reduced dy- 
namical system is described by the Hamiltonian, Hq and 
Af is taken to be the number constant. In the second 
case, Af acts as the Hamiltonian. Since, the analysis in- 
volves two infinite sums, both sums are ignored in both 
cases in the leading term analysis. In both cases, for 
unstaggered stationary localized states two permissible 
values, namely f3s and /3; of the width parameter, f3, are 
found. It is further found that for two real roots to exist 
we need v < t'criticai if v > 0. Furthermore, z^criticai is 
found to be a monotonically decreasing function of the 
parameter, a, in the first case and a in the second case. 
These parameters are defined in the text and are positive 
semidefinite. It is successfully argued from our numeri- 
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cal analysis that SLSs, characterized by the smaller width 
parameter, (3^ are stable and states characterized by /3/ 
are unstable. However, for staggered SLSs both proce- 
dures have yielded a single value, Ps and our numerical 
results indicate that these are stable localized modes of 
the system. Though both procedures yield qualitatively 
same results for the parameters of SLSs, no quantitative 
comparison is attempted here. Again, it is found that 
the problem can be exactly solved in the second case. In 
our exact solution no unstable SLS is obtained. So, the 
occurrence of unstable SLSs in this system in the leading 
term analysis should be attributed to the truncation er- 
ror. In the context of SLS in one dimensional nonlinear 
systems, this is indeed an important result. 

The formation of unstaggered and staggered SLSs are 
investigated here for v >Q. For the null value of the 
present variational procedure correctly produces SLSs of 
AL equation. Furthermore, when — > 0, it is found 
that /3s o? in the first case for both unstaggered and 
staggered SLSs. Consequently, AL stationary soliton is 
recovered in this asymptotic limit. In the second case, 
a — !■ sinh fig asymptotically as — > in the leading term 
analysis. This result is true for both unstaggered and 
staggered SLSs. The same asymptotic results are found 
in the exact analysis too, except that /3s is replaced by 
/3root- So, in the second case too the AL stationary soliton 
is the asymptotic result for v 0. Analytically also 
the same asymptotic result is obtained here. Our both 
analytical and numerical results are expected on physical 
consideration. 

In the other asymptotic analysis, a — > in the first 
case and a ^ in the second case. Again, in the second 
case there are two scenarios, the leading term analysis 
and the exact calculation. For all cases and for both 
unstaggered and staggered localized states, the known 
asymptotic form of 'J for the stable SLS are obtained 
numerically from the present variational analysis. This 
is, therefore a very important contribution of the present 
work. For unstaggered SLSs, it is found that the width of 
the state decreases with increasing v. On the other hand, 
for staggered SLSs, the width increases with increasing h' 
and vanishes a,s ly ^ 1. These results are consistent with 
the physics of the problem and reasons are given in the 
text. 

Another important aspect is the dependence of the 
width of SLSs on the position of the maximum amplitude, 
denoted by xq- It is proved in the text that Xq = or ±i. 
It is observed in our analysis that for unstaggered SLSs 
the on-site peaked SLS {xo = 0) has smaller width than 
the inter-site peaked {xq — ±i) SLS for a given value 
f > and a > 0. Our analysis also shows that for a 
given 1/ and a, the on-site peaked unstaggered SLS is en- 
ergetically more stable than the corresponding inter-site 
peaked SLS. These results are physically realistic and are 



successfully explained using the effective mass picture. It 
is found in our analysis that the existence of the P like 
mode and the ST like mode is a fundamental property 
of the system, described by IN-DNLS. It is also shown 
in this numeroanalytical method that the P like mode is 
energetically more stable than the corresponding ST like 
mode. These results constitute a very important contri- 
bution of the present work. 

It is definitely important to find exact eigenvalues and 
eigenvectors of the problem. Present analysis may turn 
out to be a useful guide for the exact calculation. In this 
analysis we have not proved that the lowest eigenvalue is 
obtained. Furthermore, the system may have more than 
one SLS type of ILM. These questions need to be properly 
investigated. Presence of impurity in the nonintegrabil- 
ity parameter, i/ may produce more stationary localized 
states and these states may interact through further ex- 
ternal perturbation. A study of this type is important in 
the transport in nonlinear systemsJ^^^J. In our calcu- 
lation, we find both P like mode and ST like mode. It 
will be interesting to know the asymptotic form of these 
modes in this model. Furthermore, the behavior of these 
modes under external perturbation should also be inves- 
tigated. Finally, it is important to find more physical as 
well as biological systems, where IN-DNLS can be used 
to study transport properties. A good candidate in this 
regard is the transport across biological membranes of 
protons through proton-wires. 



APPENDIX A: THE FORMULATION OF 
DISCRETE VARIATIONAL APPROACH FOR 
IN-DNLS 

We are dealing with a nonlinear eigenvalue problem in 
our aim to find stationary localized states of IN-DNLS 
equation, Eqs.(2.2) and (2.5) in the text. For this pur- 
pose we are employing variational formulation^S, 53, 58|. 
To implement the variational approach for this problem, 
we require the proper functional, F whose constrained 
variation will lead to Eq.(2.5)|53- We, of course know 
a constant of motion and the Hamiltonian, A/", and H 
respectively for the problem at handj^. 

Inasmuch as we know J\f, and H, using the analo- 
gous variational approach of finding eigenvalues in stan- 
dard Sturm-Liouville problems [53|. we set up the func- 
tional, F = H — A A/" where A is the Lagrange 
multiplier [H^. We then have for the variation of F, 
SF — SH — A 6J\f. For the calculation of the vari- 
ation, we transform '^n ~^ + S'^n, n G Z in the 
expression of H and Af (Eqs.(2.6) and (2.7) respectively) 
to obtain 



20 



S F 



5 H 



-2E 



- ASAf 
A + 



1+ *2 



2 



A 



(Al) 



Since, {(5^n} are arbitrary, (5F = implies that 

A (l + 'i'Di^n+i + *„-i) + 2 + A = 0. (A2) 

We note that (A. 2) is identical to Eq.(2.5) when A = 
Lu. From further analysis, we find that w is given by 
Eq.(2.12) in the text. 

For the case, where Hq = constant, the corresponding 
functional, F should be given by F = A2 Ho + 2 A/", 
where (A2 — 1) is the Lagrange multiplier. The same 
procedure will yield Eq.(2.5) if A2 



2u 



2v~ 



So, it is important to note that we can devise the re- 



quired functional to determine the eigenvalues of SLSs 
by variational approach without the formal knowledge of 
the Lagrangian. 



APPENDIX B: CALCULATION OF THE 
FUNCTION, AA(*,/3,a;o), EQ.(2.25) 

Before we proceed in this section, we cite some results 
required for the calculation|6l|. 



From (Bl) we get 



cos TT sini^ arcsinhv^l 

/(s, *,/?)) = / dy ^ — = 5-. (Bl) 

^ ' ^ coshy + (1 + 2*) 2 (1 + *) sinh2^ ^ ' 



^, ^ ^, arcsinh\/¥ d (arcsinh\/^)^ 
lim /(s,^,/3)= ^ = = -4 B2) 



In our calculation, we are using the following ansatz. 

1 



cosh/3(n — a;o) 



, n e Z. (B3) 



This ansatz has also been used in the previous 
analysis [3^. For on-site peaked and ST like localized 
states, xq = 0, and for inter-site peaked and P like 
states, xo = ±5'2l|4l|4l|43. We further write 
$2 = vj/. The function, Af{^, f3, xq) is given by Eq.(2.6) 
in the text. Now introducing (B3) in Eq.(2.6) and then 
taking partial derivative with respect to ^E", we get 



E 



n —~oo 



1 



cosh - Xo) + * 



(B4) 



We use next the famous Poisson's sum formula, Eq.(2.17) 
in the text in (B4)^3|- Thereafter, some simple algebraic 
manipulations are done to obtain 

3J\f 2 A ^ 

= -I{0,^,(3) + - ^cos(2msxo) /(s,*,/?). 

^ ^ s =1 

(B5) 



5* 



We note that (B5) is identical to Eq.(2.20) in the 
text. Furthermore, we have from the definition that 
A/'(0,/?,a;o) = 0. See Eq.(2.6) in the text. After in- 
tegration of (B5) over v?' € (0, \1/) we get Eq.(2.25) in 
the text. When = sinh^ n e Z in Eq.(2.25), 
we get A/'(5', Xo) = 2t?(5. See Eq.(3.1) in this con- 
text. This particular result has been obtained by another 
route in the literature jlj, 36]. We consider now = or 
the AL equation. Then from Eqs.(2.13), (2.20), (2.23), 
(2.24) and (2.26) we get for n e iV 



tanh n/3 coth /3 — n = 0. 



(B6) 



We see then that if n > 1, (B6) has no real nonzero /3 as 
a solution. 



APPENDIX C: APPLICATION OF THE 
METHOD OF SUCCESSIVE SUBSTITUTIONS IN 
THE INVESTIGATION OF THE FORMATION OF 
STATIONARY LOCALIZED STATES IN IN-DNLS 

Following the text, we take A/'(^', /3, xo) — 2 = 
Constant. We further write Eq.(2.25) in the form "if = 
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F{'$) where the function F{'^) is defined as 



(CI) 



and in (CI) 



/2(/3,2;o,*) = 2 /3cos[27ra;o; 



sin^(^ arcsinh\/^) 



sinh ^ 



It should be noted that only the first term in the sum 
in Eq.(2.25) is retained to obtain (CI). To obtain roots 
of the equation, 'i' = F{'^), we can use the method of 
successive substitutions. In this method at the k-th itera- 
tion, we write ^'fe+i = F('^k) with the assumption that 
limfc ^ oo ^'fe — > ^root- However, the necessary con- 
dition for this to happen is that |i^'(4'r-oot)| < l|60|. 
To explain the use of this method in the calculation, 



we shall restrict ourselves only to the first iteration with 
= sinh^ a^/p. This in turn implies 



* ^ ^-1 = sinh2[v/a2 (3 - /2(/3,xo,^'o)]. (C2) 

Note that in the calculation of roots (Eq.(3.11) in the 
text) we have taken /2 = 0. By approximating 'J by ^'i 
and furthermore keeping only the first term in the sum in 
the definition of /i(/3, t^. A, xq) (Eq. (2.21)) we get from 
Eq.(2.22) 



Ho = -4 A 4 A i^ — [1 -I- 2cos27ra;o — ^1. 

sinh/3 f3 ^ sinhf ^ 

JC3) 

We then get the permissible values of /? by setting = 



0. 
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